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Abstract 

We study the direct CP violation effect in the decay mode B + —* (j)K + . This 
decay mode is dominated by the loop induced b — > sss penguin diagram with a tiny 
contribution from the annihilation diagram. Therefore, the standard model expectation 
of direct CP violation is negligibly small. Using QCD factorization approach we find 
the CP asymmetry in the standard model to be at percent level. We consider then two 
scenarios beyond the standard model, the model with an extra vector dike down quark 
(VLDQ) and the R-parity violating supersymmetric model (RPV) and show that the 
direct CP violating asymmetry in B + — > <pK + could be as large as ~ 85% (70%) in 
VLDQ (RPV) model. 
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1 Introduction 



The study of B physics and CP violation is now at the center stage of high energy physics 
research with dedicated B factories (BABAR and Belle) already having an accumulation of 
huge data in the B sector. With these two B factories in full operation and hadronic machines 
coming up, the flavor sector of the standard model (SM) will be subjected to stringent tests 
in the near future. The major goal of these B factories is not only to test the predictions of 
the SM but also to reveal the presence of new physics (NP), if any. CP violation in B system 
has been confirmed recently at both the B factories with the measurement of sin 2(3 (where 
(3 is one of the angles of CKM unitarity triangle) in the B° — > ipKs decay [1]. It should be 
emphasized here that the current measured value of (sin 2(3)^k s by both the B factories are 
very close to each other [2]. The present world avarage on (sin2/3)^ s [2] is 

sin(2/3)^ s = 0.736 ± 0.049 , (1) 

which is in very good agreement with the SM expectation. However, this result does not 
exclude intersting CP violating new physics effects in other B decays. The rare decay 
mode B° — * (pK s , which is a pure penguin process, involving the quark level transition 
b — > sss, is one of the channels which provides powerful testing ground for new physics. The 
reason is very simple because in the SM the direct decay amplitudes for B° — > ipKs and 
B° — > <f)Ks modes have vanishing weak phase (in the Wolfenstein parametrization). Thus, 
the time dependent mixing-induced CP asymmetry in both decays is due to the weak phase 
in B° — B° mixing and expected to give the same value, i.e., sin 2(3. In the SM, the difference 
between these two measurements is expected to be very small [3] 

| sin(2/3)^ fl - sin(2/?)^ s | £ 0{\ 2 ) , (2) 

with A ~ 0.22. However, the recent measurements of sin 2(3 in the B° — > <pKs decay at 
BABAR and Belle give the average value [2] 

sm{2(3)^ Ks = -0.15 ± 0.33 , (3) 

which is about 2.7a deviation away from the corresponding measurement of (sin2/3)^ s . 
The discrepancy between the measured values of (sin 2(3)^k s and (sin 2(3)^3 may imply the 
possible existence of new physics in one of the decay modes. 

If there is new physics, then it could affect the B° — B° mixing as well as the decay 
amplitudes. But the NP contribution arising from mixing, will be same (i.e., universal) in 
both the processes . So when we compare the measurements of sin 2(3 in both the two decays 
then the NP contribution, if any in mixing, will not appear in the comparison, since it will 
affect both the decays with the same amount. Whereas, the NP contribution present in 
the decay amplitudes are nonuniversal and process dependent and can vary from process to 
process. The striking difference between these two decays is that B° — > ipKs is a tree level 
{b — > ccs) process whereas B° — > 4>Ks is a purely loop induced penguin process. So the 
NP contribution to B° — > i)K s is expected to be suppressed and it could be significant for 



2 



the loop induced process B° — > 4>K S . So, the discrepancy between the measured values of 
{sm2f3)^K s an d (sin 2{3) ( j )Ks may indicate the possibility of NP effect in the decay amplitude 
oiB°^<j)K s . 

Recently, various new physics scenarios have been explored to explain the above discrep- 
ancy [4, 5, 6, 7]. If new physics effects indeed are present in the decay mode B° — > <pK s , 
then one can expect to observe similar effects in other modes having the same internal quark 
structure. It is therefore important to explore other signals of new physics in order to corrob- 
orate this result. One way to search for new physics effects is to look for direct CP violation 
in the decay modes which are having a single decay amplitude in the standard model and 
hence expected to give zero direct CP asymmetry. An observation of nonzero direct CP 
violation in such modes is an unambiguous signal of new physics. 

At present, it appears from the current trend of data that new physics effects might be 
present in the B° — > (f)Ks mode. It is therefore interesting to see the effects of NP in the 
penguin dominated decay mode B^ — > ipK^ also (which is having the same underlying quark 
dynamics) alongwith a tiny annihilation contribution, which is our prime objective in this 
paper. The branching ratio and CP asymmetry measurements for this mode have recently 
been reported by both the Belle [8], BABAR [9] and CLEO [10] collaborations, which are 
given below as 

(9.4 ± 1.1 ±0.7) x 1(T 6 Belle, 
0.01 ±0.12 ±0.05 Belle, 
(lO.Olg;! ± °- 5 ) x 10_6 BABAR , 

0.04 ±0.09 ±0.01 BABAR, 
(5.5^-1 ± Q ^ x 1Q -6 CLE0 _ 



<pK + ) = (9.2 ±0.7) x 10" 6 , 

<pK + ) = 0.03 ±0.07. (5) 

At this point we ask the question whether the possibility of NP effect in the decay mode 
B + — > <pK + has already been ruled out ? In the rest of the paper our objective will be to 
closely examine this question and try to obtain a meaningful answer, if any. If one looks at 
the data on the direct CP asymmetry, then certainly nothing can be concluded at present. 
That is the central value is higher than the SM expectation, but as error bars are also 
significantly large, one simply cannot conclude/exclude the possibility of NP effect. So the 
present status is that, it is premature to say anything and to be more precise, on which side 
(i.e., within or outside) of the SM they are, as far as direct CP violation in B ± — > §K ± is 
concerned. To summarize, although the CP violating asymmetry is not very much larger 
than the SM expected value, but due to the presence of large error bars, decisive conclusion 
regarding the presence/absence of NP effect cannot be inferred. 

This in turn gives us enough room to explore NP effect. In fact if in future the data 
stabilize with even a few percent of direct CP asymmetry, then it may be very hard to 
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explain the same under the framework of the SM and eventually that may lead to the 
establishment of NP in this mode. Keeping this in mind we now study carefully to find the 
answer to our question that we raised earlier. 

In this paper, we intend to study the direct CP violation effects in the decay mode — > 
<f)K ± . We use the QCD factorization method [11, 12] to evaluate the relevant branching 
ratio and the direct CP asymmetry vAcp, in the SM. This decay mode B ± — > <f)K ± has 
recently been studied in the SM using QCD factorization [13, 14, 15]. However, the predicted 
branching ratio found to be well below the present experimental value. Next, we consider 
two beyond the standard model scenarios, the so called R-parity violating supersymmetric 
model and the model with an extra vector-like down quark. These two models were explored 
recently to explain the observed discrepancy in B° — > (j)K s mode [5, 6, 7]. 

The paper is organized as follows. Section II includes a general description of CP violating 
parameter in B ± — > ipK^ decays, while in Sect. Ill we analyze the decay mode B ± — > <j)K ± 
in the SM using the QCD factorization. The new physics effects from the VLDQ model 
and the RPV model are considered in sections IV and V respectively and in Section VI we 
present some concluding remarks. 



2 CP violating Asymmetry 



Here we briefly present the basic and well known formula for the CP violating asymmetry 
parameter. For charged B ± — > <pK ± decays the CP violating rate asymmetries in the partial 
rates are defined as follows : 

T(B- -> - T(B+ -> (f>K+) 

CP T(B~ -> 4>K~) + T(B+ -> <j)K+) ' 1 ' 

Without loss of generality we can write the decay amplitudes as 

A(B-^cf)K-) = X u \A u \e tS - + X c \A c \e^ , 

A(B + ^<j)K + ) = X* u \A u \e iS " + X* c \A c \e iS * , (7) 

where X q = V q bV* s , A u and A c denote the contributions arising from the current operators 
proportional to the product of CKM matrix elements \ u and A c respectively. The corre- 
sponding strong phases are denoted by 5 U and 5 C . Thus the CP violating asymmetry is given 
as 

. -2lm(X u X* c )lm(A u A* c ) 



\X U A U \ 2 + \X C A C \ 2 + 2Re(X u X*)Re(A u A*) 
2 sin 7 sin (5 U — 5 C ) 

I ^u-A u I i I \ C A C 



2 cos 7 cos(<5„ — 5 C 



(8) 



where —7 is the weak phase of A M , and A c is real in the Wolfenstein parametrization. Thus 
to obtain significant direct CP asymmetry, one requires the two interfering amplitudes to be 
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of same order and their relative strong phase should be significantly large (i.e., close to n/2). 
However, in the SM, the ratio of the CKM matrix elements of the two terms in Eq. (7) can 
be given (in the Wolfenstein parametrization) as |A U /A C | ~ A 2 V p 2 + i] 2 — 2%. Therefore, 
the first amplitude will be highly suppressed with respect to the second unless A u » A c . 
Therefore, the naive SM expectation on Acp is that it is negligibly small. This in turn makes 
the mode interesting to look for the NP in terms of large direct CP asymmetry. 

In the presence of new physics the amplitude can be written as 

A(B~ -> <j>K~) = Asm + A NP = A SM [l + r NP e<*» p ] , (9) 

where r^p = \A^p/Asm\, {Asm and A^p correspond to the SM and NP contributions to 
the B~ — > (j)K~ decay amplitude respectively) and (p^p = axg(A NP / A S m), which contains 
both strong and weak phase components. The branching ratio for the B~ — > (j)K~ decay 
process can be given as 

Br(£- -^<f>K~) = Br SM (l +r 2 NP + 2r NP cos (f) NP ) , (10) 

where Br SM represents the corresponding standard model value. 

To find out the CP asymmetry, it is necessary to represent explicitly the strong and weak 
phases of the SM as well as NP amplitudes. Although, it is expected that the SM amplitude 
\ U A U is highly suppressed with respect to its \ C A C counterpart, for completeness we will 
keep this term for the evaluation of Acp- We denote the NP contribution to the decay 
amplitude as A^p = \A^p\e l ^ n+en \ where 5 n and 9 n denote the strong and weak phases 
of the NP amplitude respectively. Thus, in the presence of NP, we can explicitly write the 
decay amplitude for B~ — > <fiK~ mode as 

A(B~ -> cf>K-) = \ u \A u \e iS " + \ c \A c \e i5 c + \A NP \e l ^ +e ^ . (11) 

The amplitude for B + — > <f)K + mode is obtained by changing the sign of the weak phases of 
the amplitude (11). Thus, the CP asymmetry parameter (6) is given as 



2iA\ u \ c A u A c \ sin7sin5 uc + \\ U A U A NP \ sin(7 + 6 n ) sm.5 un + \X C A C A NP \ sin6> n sin5 c 

\A\ 2 + 2||A U A M | ^|A C A C | cos7cos5 uc + \A NP \ cos(7 + 6 n ) cos5 Mn ^j + \\ C A C A NP \ cos 6 n cos5 cn | 

(12) 

where \A\ 2 = (\\ U A U \ 2 + |A C A C | 2 + |Aatp| 2 ) and = Si — Sj are the relative strong phases 
between different amplitudes. If we neglect the \ U A U component in the decay amplitude 
(11), the CP asymmetry will be reduced to 

2 ( | \ C A C A NP | sin 6 n sin(5 c - S n ) ) 

-4cp = ^ ? L v - (13) 

\\ C A C \ 2 + |Aatp| 2 + 2l\X c A c A NP \ cos6» n cos(5 c - S n )\ 

Having obtained the direct CP asymmetry parameter Acv in the presence of new physics 
we now proceed to evaluate the same within and beyond the SM, in the following sections. 
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3 CP Violation in B ± 0i^ ± process in the SM 

To study the CP violation effects in B ± — > <pK ± process, first we present the SM amplitude 
and find out the branching ratio and CP asymmetry parameter. In the SM, the decay 
process B ± — > (pK^ receives dominant contribution from the quark level transition b — » sss, 
which is induced by the QCD, electroweak and magnetic penguins and a tiny annihilation 
contribution. The effective AB = AS = 1 Hamiltonian relevant for the process under 
consideration is given as [11] 



H eff 



G f 

V2 



li 



r- 



(14) 



where q = u, c. Oi^, O3, • • • , 0% and O7, • • • , O10 are the standard model tree, QCD and 
EW penguin operators, respectively, and On is the gluonic magnetic operator. The values 
of the Wilson coefficients at the scale /i m b in the NDR scheme are given in Ref. [16] as 

Ci = 1.082, C 2 = -0.185, C 3 = 0.014, C 4 = -0.035 , 

C 5 = 0.009 , C 6 = -0.041 , C 7 = -0.002a , C 8 = 0.054a , 

C 9 = -1.292a , C10 = 0.263a , C n = 0.052 . (15) 

We use QCD factorization [11, 12] to evaluate the hadronic matrix elements. In this 
method, the decay amplitude can be represented in the form 



{</>K-\Oi\B-) = (</>K-\Oi\B- 



'fact 



(16) 



where (0K~|Oj|-B~)f act denotes the naive factorization result and Aq CD ~ 225 MeV, the 
strong interaction scale. The second and third terms in the square bracket represent higher 
order a s and Aqcd/^ corrections to hadronic matrix elements. We shall closely follow the 
approaches [13, 14, 15, 17] to reanalyze our QCD factorization calculation. In a recent paper 
Beneke and Neubert have presented the latest QCD factorization calculation [18], where 
the authors have considered 4 different schemes to find a suitable set of parameters which 
can explain the observed data. Furthermore, it has been shown in that paper that the last 
scheme, i.e., S4 (scheme 4) explains the observed data on branching ratios and CP violation 
parameters to a good accuracy. Of course it should be mentioned here that even the scheme 
S4 also fails to explain the problematic penguin dominated B — > <pK modes [18, 19]. 

In this analysis, however, we will restrict ourselves to the default QCD factorization 
(by default we mean the QCD factorization result without any parameter readjustment (for 
more details the reader is referred to Ref. [18]). This is done so, because of the fact that 
one should not just ignore the possibility of NP effect. Furthermore, we do not see why 
the parameter readjustment is the only possibility to match the data (eventhough it may 
be a logistic option and in the long run may even come out to be true). So, in order to 
visualize the NP effect, we confine ourselves to the default values of the QCD factorization 
and allow NP to fill up the gap between the experimental data and the SM values rather 
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than readjusting the parameters to find a good agreement. In doing so, we can accommodate 
the beyond the SM scenarios, if any, to explain the possible discrepancy. Nevertheless, we 
will compare the results of Ref. [18] wherever necessary to that of our work. It should be 
emphasized here that our speculation regarding the presence of NP may not be unrealistic. 
As a matter of fact the present data in this sector do not seem to comply with the SM 
expectations, already mentioned earlier, as in the case of B° — > 4>Ks decay. 

In the heavy quark limit the decay amplitude for the B~ — > (f>K~ process, arising from 
the penguin diagrams, is given as 



a p {b- - 4>k~) = % £ v qb v; 

q=u, c 



a q 3 + a q 4 + a q 5 - - (a q 7 + a q 9 + a q 10 ) + a q 10a 



X 



(17) 



where X is the factorized matrix element. Using the form factors and decay constants defined 
as [20] 



(K-(p K )\*fb\B-(p B )) 



(Pb+PkT- 



m B - m K » 



Fi{q 2 ) 



+ 



m R — m 



Mq,e)\*f8\0) 



(18) 



we obtain 



X = (K-( Px )|s 7m (1- 75 )&|S-( Pb ))(0( ? , e )|^(i_ 75 ) s |o) 
= 2F 1 B ^ K (ml)Um^(e-p B ). 



(19) 



The coefficients a^'s which contain next to leading order (NLO) and hard scattering correc- 
tions are given as [13, 17] 



C A a s C F 
° 3+ iV+^A^' 



n C 3 a s C F 
° 4+ iV AT 



CslF^ + G^+G^st,)] , 



+ C.G^s,) + (C 4 + C 6 ) J2 G^sf) + CuG^u 

f=u 

= «8 = C. + £ + £%CW-F,-12), 



no 



a% = a% 



On — CLq 



no 



N Air N 
C 9 a s C F 
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n q 
a lOa 



a s C F 



/ f=u 



+ Cg- (esG^Ss) + e 6 G^(s 6 )) 



(20) 



where q takes the values u and c, iV = 3, is the number of colors, Cf = (N 2 — 1)/2JV. The 
internal quark mass in the penguin diagrams enters as s/ = nij/m 2 . The other parameters 
in (20) are given as 



GJs) 



-12^^-18 + 4 + ^, 

nib 

dx g{x)(t) (t> (x) , 
X \nx — Zm , 



1 — x 

47T 2 f K f. 



B 



1 dz 



N F^ K (0)m 2 B Jo z 



4> B {z) f —<p K (x) f —4><t>{y) , 

Jo x Jo v 



2 4 , A* „ 
- ln-^ + 4 

3 3 nib 



dx 



x) 



X 



Ga 



,11 



/ c£u it(l — u) In [s — u(l —u)(l — x) 
Jo 

r 1 2 

- / dx- 4>Ax) . 

Jo 1 — x 



(21) 



The light cone distribution amplitudes (LCDA's) at twist two order are given as 

4>b(x) = N B x 2 (l — :r) 2 exp^ — 



2 2 

m B x 



<t>K,4>{ x ) = M 1 ~ x ) > 



(22) 



where iV^ is the normalization factor satisfying Jq 1 dx<pB{x) = 1 and = 0.4 GeV. The 
quark masses appear in G(s) are pole masses and we have used the following values (in 
GeV) in our analysis; 



m>u = nid = m s = 0, m c = 1.4 mj = 4.8. 
The contributions arising from the annihilation diagrams are given [14, 15] as 



G 



—f=fBfKf<t> 2 

V2 m% 



VubV* s {b 2 
+ V cb V: s {b^K-) + bT^K-)} 



K~ 



\K- 



(23) 



The parameters bi [15] are as 



br(<f>,K~) 



N 2 
N 2 



C 3 4(0, K~) + C 5 4(0, K~) + (C 5 + iVC 6 ) 4(0, if~ 
C 9 Ai(0, A-) + C 7 4(0, if") + (C 7 + iVC 8 ) 4(0, if)" 



(24) 



where 



4(0, if") 
4(0, if-) 
4(0, 



7Ttt s 







/ dx 


^(x) / dy X (y) X 


/o 


JO 



+ 



= 7ra 



vra s jf da; ^(x) dy <f> K (y) x r x^y 
dr (x) jf x (y) > 



2(1 -y) 



y(l - x)(l - x(l - y)) ' 
2(2 -x) 
y(l-a:) 2 ' 



(25) 



These integrals contain the divergent end point integrals. Assuming SU(3) flavor symmetry 
and symmetric LCDA's (under x <-> (1 — x)), the weak annihilation contributions can be 
parametrized as 



4(0, if~) - 187ra s 



7T 



X A -4+ — 



4(0, K~) ~ 7r« s r x [2tt 2 + 6(X^ - 2X A ) 



2X A — X A 



(26) 



where = Jq dx/x parameterizes the divergent end point integral and r x = 2m 2 K /m b {m s + 
m u ) denotes the chiral enhancement factor. It should be noted that the quark masses in the 
chiral enhancement factor are running quark masses and we have used their values at the b 
quark mass scale as m b {m b )=AA GeV, m s (m b )=90 MeV and m u {m h ) = 3.2 MeV. Thus we 
obtain the total amplitude as (in units of 1CT 4 ) 



A(B~ -> 0if-) 



Gf 
'V2 



(e-Ps) 



V ub V* s (45.022 + i 47.788) + V cb V c * s (51.407 + i 28.054) 



G F 

~T2 {e ' PB) 



V ub V: s (65.656 e m °) + V cb V c * s (58.564 e m °) 
The branching ratio can be obtained using the formula 



(27) 



0if" 



8nm 2 



\A{B~ ^<f>K-)/{e-p B ) 



(28) 



where p cm is the momentum of the outgoing particles in the B meson rest frame. 

We have used the following input parameters. The value of the form factor at zero recoil 
is taken as F^ K (0) = 0.38, and its value at q 2 = m 2 , can be obtained using simple pole 



dominance ansatz 



as F{ 



B-*K, 



0.39. The values of the decay constants are as 
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0.233 GeV, f B = 0.19 GeV, iW=0.16 GeV, and the lifetime of B~ meson t b - = 1.674 x 10~ 12 
sec [21]. For the CKM matrix elements, we have used the Wolfenstein parameterization and 
have taken the values of the parameters A = 0.819, A = 0.2237, p = 0.224 and rj = 0.324. 
With these input parameters, we obtain the branching ratio in the SM as 

Br SM (B- -> <j)K-) = 5.99 x 10~ 6 , (29) 

which lies quite below the present experimental limit (5). The corresponding CP asymmetry 
parameter Aqp is found to be 

A C p = 0.011 , (30) 

which is also below the present central experimental value (5). However, our predicted 
branching ratio is in close agreement with the various earlier calculations [14, 15] and with 
that of the default value in [18]. The small difference between the predicted values in these 
analyses is due to difference in the used input parameters (i.e quark masses, decay constants 
and CKM parameters etc.). It should be noted here that the obtained value of branching 
ratio is slightly above the predicted value of Ref. [13], where they have not considered the 
annihilation contributions. The Aqp value obtained here is in agreement with the previous 
calculation of Beneke and Neubert [18]. It should be pointed out here that even the best 
scheme S4 [18], which predicts close to the branching ratio obtained by the experiments, 
predicts even smaller Aqp (0.007). Now if we consider the error bars in the data, certainly 
upto 10 % of Aqp is allowed by the present data. Such a large Aqp, if established later, 
certainly cannot be accommodated in the SM. Thus in turn this brings us to the door of NP 
to provide some meaningful explanation. It can be seen from Eq. (27) that since \A U \ ~ \A C \ 
and also their strong phases (S u and S c ) are nearly equal, the observed CP asymmetry is quite 
small. Furthermore, as discussed in section II, since \A U \ ~ \A C \ and V u bV* s /V c bV* s ~ 2%, we 
will neglect the \ U A U term in the decay amplitude (11) when we consider the contributions 
from beyond standard model scenarios. 



4 Contribution from the VLDQ Model 

Now we consider the model with an additional vector like down quark [22]. It is a simple 
model beyond the SM with an enlarged matter sector with an additional vector-like down 
quark D4. The most interesting effects in this model concern CP asymmetries in neutral B 
decays into final CP eigenstates. At a more phenomenological level, models with isosinglet 
quarks provide the simplest self consistent framework to study deviations of 3 x 3 unitarity 
of the CKM matrix as well as allow flavor changing neutral currents at the tree level. The 
presence of an additional down quark implies a 4 x 4 matrix V ia (i = u, c, t, 4, a — d, s, b, b'), 
diagonalizing the down quark mass matrix. For our purpose, the relevant information for 
the low energy physics is encoded in the extended mixing matrix. The charged currents 
are unchanged except that the Vckm is now the 3x4 upper submatix of V . However, the 
distinctive feature of this model is that FCNC enters neutral current Lagrangian of the left 



10 



handed downquarks : 

9 



2 cos 9 W 
with 



u L il^u Li - d La U aP 'fd L p - 2 sin 2 6 W J» m , (31) 



U*fi= E = 5 a(3 - V; a V, p , (32) 

i=u,c,t 

where U is the neutral current mixing matrix for the down sector which is given above. As 
V is not unitary, U ^ 1. In particular its non-diagonal elements do not vanish : 

£4/3 = -OW0 fora^/3, (33) 

Since the various U a p are non vanishing they would signal new physics and the presence 
of FCNC at tree level, this can substantially modify the predictions for CP asymmetries. 
The new element U sb which is relevant to our study is given as 

U sb = V: s V ub + V c * s V cb + V*V tb . (34) 

The decay modes B ± — > <pK ± receive the new contributions both from color allowed and 
color suppressed Z-mediated FCNC transitions. The new additional operators are given as 

Z-FCNC = _ lb)ha][ s pl ^C S y - C S Al ,) S p] , 

Ot FCNC = [8ff(l - 7fc)&«][*«7/.(Cv - C s Al s)sp\ , (35) 

where C v and C S A are the vector and axial vector Zss couplings. Using Fierz transformation 
and the identity (C v - C% lb ) = [{C s v + C S A ){1 - 7s ) + (C v - C S A )(1 + 75 )]/2, as done in [7], 
the amplitude, including the nonfactorizable contributions, is given as 

A V ^{B~ - 4>K~) = ^U sb 2X (c v + f ) (l + ^~F,) . (36) 

The values for C v and C% are taken as 

C^ = ~ + ^ 2 9 w , C' A = -\. (37) 

Now using sin 2 6^=0.23, alongwith U sb = \U sb \e l9n , where 9 n is the new weak phase of VLDQ 
model and 

\U sb \ ~ 1 x 1(T 3 , (38) 

we find the amplitude 

A VLDQ (B- -> cf)K-) = -^=e Wn (e • p B ) (1.683 - i 0.140) x 10" 4 . (39) 

v2 

Thus we obtain the strong phase in VLDQ model to be ~ —5° and r NP as 

r NP = 0.7 . (40) 
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From Eq. (40), it should be noted that the NP amplitude is of the same order as the SM 
amplitude and hence large CP asymmetry is expected in this model. Now plotting the 
branching ratio (9) versus 4>np (the red line in Figure-1), we see that the observed data can 
be easily accommodated in VLDQ model. The direct CP asymmetry ^4cp (Eq. (13) vs. 6 n 
is plotted in Figure-2, for 5 C = 29° and 5 n = —5°, and it is also (the red line in Figure-2) in 
compatible with the present experimental data. Furthermore, it should be noted from the 
Figure-2 that VLDQ model can accommodate the CP violating asymmetry upto 85%. 




50 100 150 200 250 300 350 

0NP 



Figure 1: Branching ratio of B~ — > <pK~ process (in units of 10~ 6 ) versus the phase (p^p (in 
degree). The red and blue curve represent the results of VLDQ and RPV model respectively. 
The horizontal solid line is the central experimental value whereas the dashed horizontal lines 
denote the error limits. 



5 Contribution from R-Parity violating supersymmet- 
ric model 

We now analyze the decay mode in the minimal supersymmetric model with R-parity vio- 
lation. In the supersymmetric models there may be interactions which violate the baryon 
number B and the lepton number L generically. The simultaneous presence of both L and 
B number violating operators induce rapid proton decay, which may contradict strict exper- 
imental bound. In order to keep the proton lifetime within experimental limit, one needs to 
impose additional symmetry beyond the SM gauge symmetry to force the unwanted baryon 
and lepton number violating interactions to vanish. In most cases this has been done by 
imposing a discrete symmetry called R-parity defined as, R p = (— l)( 3B + L+25 ') ) where S is 
the intrinsic spin of the particles. Thus the .R-parity can be used to distinguish the particle 
(Rp=+1) from its superpartner (R p = —1). This symmetry not only forbids rapid proton 
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Figure 2: The direct CP asymmetry (.Acp) m the B — > (j)K process versus the new weak 
phase n (in degree). The red and blue curve represent the results of VLDQ and RPV model. 
The horizontal solid line is the central experimental value whereas the dashed horizontal lines 
denote the error limits. 

decay, it also renders stable the lightest supersymmetric particle (LSP). However, this sym- 
metry is ad hoc in nature. There is no theoretical arguments in support of this discrete 
symmetry. Hence it is interesting to see the phenomenological consequences of the breaking 
of R-parity in such a way that either B and L number is violated, both are not simultane- 
ously violated, thus avoiding rapid proton decay. Extensive studies has been done to look 
for the direct as well as indirect evidence of R-parity violation from different processes and 
to put constraints on various R-parity violating couplings. The most general i?-parity and 
Lepton number violating super-potential is given as 

1 

2' 

where, i,j, k are generation indices, Lj and Qj are SU(2) doublet lepton and quark superfields 
and E%, D% are lepton and down type quark singlet superfields. Further, Xijk is antisymmetric 
under the interchange of the first two generation indices. Thus the relevant four fermion 
interaction induced by the R-parity and lepton number violating model is [23] 

(s«7 M (l - 7s)s/3) (%7^(1 + 7s)M 



Wy = -XijkULjEl + X' ijk LiQjDl , (41) 



+ (^£2X^(1 + 7s)s/3) - 75)& 



(42) 



It should be noted that, the factorized matrix elements of both the operators in Eq. (42) 
are same because of 

(K-\s^(l-j 5 )b\B-) = (K-\sj"(l + j 5 )b\B-) = (K-\s^b\B~) , 

- 76 ) a |0) = (01^(1 + 75 >|0) = (<j)\s^s\V>) . (43) 
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So, including the nonfactorizable effects, we obtain the R-parity violating contribution to 
the decay amplitude as 



f a s Cf , ^ 

iv + Slv<- F *- 12 ) 



(44) 



where the summation over % = 1, 2, 3 is implied. Now considering the values of R-parity 
couplings from [5] as 

where 6 n is the new weak phase of R-parity violating couplings and R ~ 1.5 x 10~ 3 for the 
mass scale M ~ m/> = 100 GeV. 

Thus we obtain the R-parity violating amplitude as 

A RPV (B -> 07^) = -(e • p B )e ie " (2.462 + % 0.512) x 10" 9 , (46) 

with the strong phase as ~ 12°. Thus, the new physics contribution r NP is found to be 

r NP = 1.28 . (47) 

Plotting the branching ratio (9) vs. 4> NP we can see that the observed branching ratio can 
be easily accommodated in the RPV model (blue curve in Fig-1). Similarly the observed 
direct CP asymmetry can also be explained by RPV model, (blue curve in Figure-2). As 
can be seen from Fig-2, the RPV model can accommodate CP asymmetry upto 70%. 



6 Conclusion 

The measured branching ratio and direct CP asymmetry parameter in the B ± — > 4>K ± mode 
are not in agreement with the SM expectations. The direct CP asymmetry in this mode in 
the SM is expected to be very small (at 1% level) but the measured value is higher than 
this. But this should be treated with caution since the error bars are still very large in 
comparison to the central value and definitive conclusion cannot be obtained. On the other 
hand there is a very wild speculation and/or a general belief that exist in the literature is that 
there might be some kind of new physics effect present, in the decay amplitude of the pure 
penguin process B° — > <f>K s , to account for the 2.7a deviation of the measured (sin2/9)^x s 
from (sin2/3)^ s . So it is quite obvious that presence of NP may not be ruled out on its 
charged counterpart process B ± — > <pK ± having the same underlying quark dynamics. 

In fact there are various ways to test the presence of new physics that exist in the litera- 
ture. One of them is to find nonzero direct CP asymmetry in modes which are dominated by 
a single decay amplitude in the SM and hence expected to give zero direct CP asymmetry. 
The mode B ± ^>K ± falls into that class, where the Aqv is expected to be negligibly small. 
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But the present data do not seem to respect the same and may reveal hidden NP effects in 
future. In order to look for NP, we first carefully reanalyzed this mode (B^ — > ifiK^) in the 
default QCD factorization approach. We also took into account of the small annihilation 
contribution in our calculation and found that the branching ratio and Acp are deviated 
from the present experimetal data. However, our results are in close agreement with the 
latest theoretical calculation [18]. 

We took this opportunity to incorporate the possible new physics scenarios to explain the 
same. In order to obtain significant direct CP asymmetry both the interfering amplitudes 
(i.e., the SM and NP amplitudes) should be of the same order and the relative phase differ- 
ences between them should be nonvanishing. We then introduce two beyond the standard 
model scenarios in turn to explain the observed branching ratio and Acv- These are the 
VLDQ model and the RPV supersymmetric model. These two models are shown in the 
literature that they can successfully explain the possible discrepancy of B° — > <pKs result 
with that of the SM expectation. We find that r NP , which is the ratio between the NP and 
SM amplitudes, to be 0.7 and 1.28 for VLDQ and RPV model respectively. As can be seen 
now from the Figures-1 and 2, that these two models can also accommodate the possible 
deviation of branching ratio and CP asymmetry parameter from the SM values for the decay 
under consideration. Eventually, it turns out that direct CP asymmetry upto ~ 85% (70%) 
can be accommodated in the VLDQ (RPV) model, if nature is so obliged to give such a 
large value. In fact more precise data in future will lead us to definitive conclusions and/or 
narrow down some NP scenarios. 

To conclude, the branching ratio and direct CP asymmetry parameter (.Acp) a s obtained 
in Belle and BABAR indicate deviation from that of the SM expectation in the B ± — > <fiK ± 
decay. If the trend remains the same, as it is now, then it may lead to the establishment of 
the presence of new physics in penguin dominated (6 — > sss) B decays. At present it is too 
early to say anything in favor of or rule out the possibility of the existence of new physics in 
this sector. At the same time keeping in mind all these we should keep every option open, 
explore various possibilities and hope that B factory data will reveal new physics in the near 
future. 
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